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Abstract. 

We used photometric data from the WASP (Wide-Angle Search for Planets) survey to explore 
the possibility of detecting eclipses and transit signals of brown dwarfs, gas giants and terrestrial 
companions in close orbit around white dwarfs. We performed extensive Monte Carlo simulations 
and we found that for Gaussian random noise WASP is sensitive to companions as small as the 
Moon orbiting a V ^12 white dwarf. For fainter stars WASP is sensitive to increasingly larger 
bodies. Our sensitivity drops in the presence of co-variant noise structure in the data, nevertheless 
Earth-size bodies remain readily detectable in relatively low S/N data. We searched for eclipses 
and transit signals in a sample of 194 white dwarfs in the WASP archive however, no evidence for 
companions was found. We used our results to place tentative upper limits to the frequency of such 
systems. While we can only place weak limits on the likely frequency of Earth-sized or smaller 
companions; brown dwarfs and gas giants (radius~ Rj up ) with periods <0.2 days must certainly be 
rare (< 10%). More stringent constraints requires significantly larger white dwarf samples, higher 
observing cadence and continuous coverage. The short duration of eclipses and transits of white 
dwarfs compared to the cadence of WASP observations appears to be one of the main factors limiting 
the detection rate in a survey optimised for planetary transits of main sequence stars. 
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INTRODUCTION 

The transit technique involves searching for periodic dips in stellar light-curves due to 
the orbital revolution of a transiting body, blocking a fraction of the stellar light. For 
a given planetary radius, the transit depth (<5) is proportional to (R p //R*) 2 . Therefore, 
planets orbiting solartype stars have extremely shallow eclipses, blocking ~1% of the 
light for a giant planet and ~0.01% of the light for an Earth-sized planet. Current 
ground-based wide-field surveys can achieve the necessary photometric accuracy of 
better than 1%, only for the brightest stars (V ~9-12 in the case of WASP), so the bulk 
of the planets discovered by transit surveys around main-sequence stars have radii in 
the range R p i ~0.8— 1.8 R; Mp . A strong advantage over main sequence star primaries is 
offered by white dwarf stars. White dwarfs (WD) are compact degenerate objects, with 
approximately the same radius as the Earth, and represent the final stage of evolution of 
main-sequence stars with masses < 8M (i.e. ~97% of all stars in our Galaxy). Any 



sub-stellar or gas giant companion in orbit around a white dwarf will completely eclipse 
it, while bodies as small as the Moon will have relatively large transit depths (~ 3%), 
with the only caveat being that it remains unclear as to whether any such systems 
survive beyond the latter stages of stellar evolution. 

Sub-stellar companions to white dwarfs are rare (Farihi et al. 2005 using 2MASS 
estimated that < 0.5% of WDs have L dwarf companions). At the time of writing only 
three wide white dwarf + brown dwarf (WD+BD) systems have been spectroscopically 
confirmed, GD 165 (Becklin and Zuckerman 1988), PHL5038 (Steele et al. 2009), 
and LSPM 1459 + 0857 AB (Day-Jones et al. 2010) and two detached, non-eclipsing, 
short-period WD+BD systems are currently known, WD0137 - 349 (Maxted et al. 2006, 
Burleigh et al. 2006, P^116mins), and GD1400 (Farihi and Christopher 2004, Dobbie 
et al. 2005, Burleigh et al. 2010, P^9.9h). The latter, is currently the lowest mass 
(~50M JM/7 ) object known to have survived CE evolution. Although infrared surveys 
such as UKIDSS, VISTA and WISE, and observatories such as Spitzer hope to reveal 
many more such binaries, they remain difficult to identify either as infra-red excesses or 
through radial velocity measurements. In addition the detection of a significant number 
of eclipsing WD+BD binary systems might help uncover the hypothesised population 
of 'old' cataclysmic variables (CVs) in which the companion has been reduced to 
a sub-stellar mass (e.g. Patterson 1998; Patterson et al. 2005; Littlefair et al. 2003). 
These systems are undetectable as X-ray sources and difficult to identify in optical 
and infra-red surveys. Littlefair et al. (2006) confirmed the first such system through 
eclipse measurements, while Littlefair et al. (2007) showed that another eclipsing CV, 
SDSS J150722.30 + 523039.8, was formed directly from a detached WD+BD binary. 

Several theoretical studies discuss post-main sequence evolution of planetary systems 
and show that planetary survival is not beyond possibility (Duncan and Lissauer 1998; 
Debes and Sigurdsson 2002; Burleigh et al. 2002; and Villaver and Livio 2007). Radial 
velocity observations of red giants indicate that planets in orbits beyond the red giant's 
envelope can survive stellar evolution to that stage (see Frink et al. 2002; Hatzes et al. 
2005, Sato et al. 2003). Moreover, Silvotti et al. (2007) reported the detection of a 
~3Mj up planet orbiting an extreme horizontal branch star, and Mullally et al. (2008) 
found convincing evidence of a 2Mj up planet in a 4.5 year orbit around a pulsating 
WD. Furthermore, Beuermann et al. (2010) reported the detection of two planetary 
companions (M c =6.9Mj up and Mj=2.2My M/ ,) in the post common envelope binary 
NN Ser (ab) via measurements of a light-travel-time effect superposed on the linear 
ephemeris of the binary; showing that planets do survive stellar evolution. 

Short-period rocky companions to white dwarfs may seem less likely. Villaver and 
Livio (2007) suggested that planets in orbit within the reach of the AGB envelope will 
either evaporate or in rare cases, more massive bodies may accrete mass and become 
close companions to the star. Planets in wide orbits that escape engulfment by the 
red giant or asymptotic giant will move outwards to conserve angular momentum (as 
described by Jeans 1924). Duncan and Lissauer (1998) found that for WD progenitors 
experiencing substantial mass loss during the AGB phase, planetary orbits become 
unstable on timescales of <10 8 year. Debes and Sigurdsson (2002) found that the mass 



loss is sufficient to destabilise planetary systems of two or more planets and that the 
most likely result is that one planet would be scattered into an inner orbit (occupied, 
before the RGB phase, by a 'now evaporated' inner planet), while the other would either 
be boosted into a larger orbit, or ejected from the system altogether. 

The above scenario provides a plausible explanation for the recent detection of silicate- 
rich dust discs around a growing number of white dwarfs at orbital radii up to ~1R 
(e.g. Reach et al. 2005; Farihi et al. 2007, 2008; Jura 2003). Jura (2003) suggests that the 
formation of dust discs around white dwarfs is most probably due to the tidal disruption 
of an asteroid or larger body which has strayed too close to the parent star. (Jura et al. 
2009) suggest that the disc around GD362 originated from the tidal destruction of a 
single massive body such as Callisto or Mars. 

The detection of short period sub-stellar and planetary mass companions to white 
dwarfs, will open an exciting chapter in the study of exoplanet evolution, constraining 
theoretical models of CE evolution and helping us to understand the ultimate fate of hot 
Jupiter systems as well as the fate of our own solar system in the post main-sequence 
phase. Here we present some of the results of our study which investigated the detection 
limits for transiting sub-stellar and terrestrial companions in close orbits around white 
dwarfs (for more details see Faedi et al. 2010). 

In §2 we discuss the characteristics of the transit signals, the parameter space investi- 
gated and our detection method. In §3 we analysed a sample of 194 WDs in the WASP 
archive. Finally in §4 we discuss our conclusions. 

CHARACTERISTICS OF THE TRANSIT SIGNAL 

A transit signal is described by its duration, its depth and its shape. Extra-solar planets 
transiting main-sequence stars show signals characterised by an ingress, a flat bottom 
and an egress, with a typical durations of 2-3 hours and depths of about 1% (see for 
example Collier Cameron et al. 2007; Simpson et al. 2010; Barros et al. 2010; and 
Faedi et al. 2010). We modelled the synthetic dataset assuming circular orbits and fixed 
stellar parameters. We considered a typical 1 Gyr old carbon-core white dwarf of mass 
M* = O.6M and radius R* = 0.013R©. We explored the detectability of planetary 
transits across the two-dimensional parameter space defined by the orbital period and 
the planet radius. Our simulations cover companions ~ 0.3R® < R p i < 12R©, and 
orbital periods in the range P ~ 2 hours to 15 days (equivalent to orbital distances 
between a ~ 0.003 and 0. 1 AU). We chose the minimum orbital period to yield an orbital 
separation close to the Roche radius of the WD, and the maximum period in order to have 
a reasonable chance of detecting five or more transits in a typical WASP season of 150 
day. 

In Figure 1 we show the probability that a given system will transit, and the duration of 
such transit across the parameter space defined above. It is evident from these diagrams 
that the signatures of transits of white dwarfs by typical planet- sized bodies will be 
rather different than those seen for typical transiting hot Jupiters. In particular the transit 
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FIGURE 1. Contours of constant transit probability (left), and duration (right) in the parameter space 
defined by orbital period and planetary radius. The transit probability is expressed in percentage values. 
The transit duration is expressed in minutes. 



duration is much shorter, from ~l-30 min for companions with sizes ranging from 
Moon-size to Jupiter-size, compared to 2-3 hours for a typical hot Jupiter. In addition, 
Figure 2, left-panel, shows that the transit depths are much larger, from around 3% for a 
Moon-sized to 100% for any companion larger than the Earth, compared to ~1% for a 
hot Jupiter. 



Synthetic WASP light-curves 

The synthetic light-curves were generated using the time sampling of a typical WASP 
survey field, and with statistical S/N representative of magnitude spanning the range of 
brightness of WDs in the WASP survey. The corresponding photometric accuracy of 
WASP over this range is ~1% to 10%. Because WASP data show residual covariant- 
noise structure we have tested the transit recovery rate in the case of both uncorrelated 
"white" noise and correlated "red" noise. To cover the orbital period-planet radius 
parameter space we selected seven trial periods spaced approximately logarithmically 
(P = 0.08, 0.22, 0.87, 1.56, 3.57, 8.30 and 14.72 days), and five planet radii R p i = 10.0, 
1.0, 0.6, 0.34 and 0.27 R®. We modelled the set of synthetic light-curves by injecting 
fake transit signals into phase-folded light-curves at each trial period with a random 
transit epoch to in the range < to < P. Because in the case of a WD host star con- 
sidered here, the ingress and egress duration is typically short compared to cadence 
of the WASP survey (8-10 minutes), we ignored the detailed shape of the ingress and 
egress phases and modelled the transit signatures as simple box-like profiles. Figure 
2, right-panel shows two examples of our simulated transit light-curves. The top panel 
shows the synthetic light-curve of an hypothetical eclipsing WD+BD binary system 
with an orbital period of P = 116 mins, similar to WD0137 — 349 (a non-eclipsing 
system, Maxted et al. 2006). The lower panel shows the simulated transit light-curve for 
a rocky body of radius 1.2R© in a 5 hr orbit. 
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FIGURE 2. Left-panel: contours of constant transit depth; right-panel: examples of synthetic light- 
curves. Top, an eclipsing BD in orbit with 2 hr period. Bottom, a 1 .2R e companion in 5 hr orbit. 



Detecting transit signals 

To recover the transit signals from the synthetic light-curves we used an implemen- 
tation of the box-least- squares (BLS) algorithm (Kovacs et al. 2002) commonly used 
to detect transits of main sequence stars. To ensure that the BLS search was sensitive 
across the expected range of transit durations, we chose to search a grid of box widths 
Wb = {1 , 2,4, 8, 16, 32} minutes, covering the range in transit durations over most of our 
parameter space (Figure 1). In addition, we used an optimised version of the BLS code 
which best accounted for the shape, duration and depth of the signals investigated in this 
work (Faedi et al. 2010). 

We used the Signal Detection Efficiency (SDE) metric defined in Kovacs et al. (2002) 
to assess the likely significance of a peak in a BLS periodogram. We evaluated SDE as 
follows: 

SDE = S P eak ~ S 

where S pea k is the height of the peak, and S and c$ are measures of the mean level 
and scatter in the noise continuum of the periodogram. A detection is represented 
by the highest peak in the BLS power spectrum. We regard as a match any trial in 
which the most significant detected period is within 1 % of being an integer fraction 
or multiple from 1 /5 x to 5 x the injected transit signal. Details of the algorithm false 
alarm probability can be found in Faedi et al. (2010). The results of our simulations 
are illustrated in Figure 2, left-panel in the case of a V ~12 magnitude WD for light- 
curves with red noise. It is evident from Figure 2; Figure 6 and Table 1, 2, and 3 from 
Faedi et al. (2010), that transiting companions are essentially undetectable at our longest 
trial periods (8.30 and 14.72 days) in a WASP-like survey; the transits are too short in 
duration and too infrequent to be adequately sampled. In addition, we found that for 
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FIGURE 3. Left-panel: recovery rate for simulated transit signals injected into synthetic light-curves 
of a white dwarf of magnitude V ~ 12; right-panel: upper-limits on companion frequency (95%) folding- 
in the detectability of transiting systems in a WASP-like survey. In both panels values are expressed in 
percent. 



idealised photon-noise-limited cases, objects as small as Mercury could be detected to 
periods of around 1 .5 d, and the Moon for periods less than 1 d. Once red noise is added, 
Moon-sized companions become almost undetectable. However, for companions around 
1R® and larger there is a good chance of detection out to periods of around 4 days. 
Our key conclusion from these simulations is that for the case of transits of white dwarfs 
the degree of photometric precision delivered by a survey is of somewhat secondary 
importance compared to a high cadence and continuous coverage. For planet-sized 
bodies individual transits will be quite deep and readily detectable in data of moderate 
photometric quality, however it is the short duration of the transits that is the main factor 
limiting the transit detection rate in surveys optimised for main sequence stars. 

SEARCHING FOR TRANSIT SIGNALS IN WASP SURVEY DATA 

Encouraged by the results of our simulations we selected a sample of 194 WDs (with 
V < 15) which have been routinely monitored by WASP through the 2004 to 2008 
observing seasons, and performed a systematic search for eclipsing and transiting sub- 
stellar and planetary companions. We selected the sample by cross-correlating the WASP 
archive with the McCook & Sion catalogue (McCook and Sion 2003). In addition to our 
authomated search, we have inspected each of the individual light-curves by eye. In 
both searches we found no evidence for any transiting and eclipsing companions. We 
have used this null result together with the results of simulations to estimate an upper- 
limit to the frequency of such close companions for the sample of WDs considered in 
this study. 

In order to estimate an upper limit to the frequency of close sub-stellar and planetary 
companions to white dwarfs, we used the detection limits derived from our simulations 
and the results obtained from the analysis of the sample of 194 white dwarfs. Although 
our complete sample numbers N = 194 stars, only a fraction p tr {R p i,P) will exhibit 



a transit, and of those only a fraction Pdet(RphP) would be detectable in a WASP- 
like survey. Both these factors act to reduce the total number of transiting companions 
detected in the survey, or in the case of a null result will tend to weaken the constraints 
that can be placed on true companion frequency by such a survey. We incorporate these 
factors and we modified our effective sample size as N' = N x p tr [R p i , P) x pj et (R p [ , P) . 
We combined the magnitude- specific p^et maps obtained from our simulations for WDs 
of magnitude V ~ 12, 13, 15 into a single map by interpolating/extrapolating according 
to the magnitude of each object in our sample and combining these to form an averaged 
map which can be folded in to our calculation of the upper-limits. The resulting limits 
corresponding to the 95% of the integrated probability, are shown in the right panel of 
Figure 3. Our results show that for rocky bodies smaller than the size of Mercury no 
useful upper limits to the frequency of companions to white dwarfs can be found, and 
that for Earth-sized companions only weak constraints can be imposed. However, it does 
suggest that objects the size of BDs or gas giants with orbital periods P< 0. 1 — 0.2 days 
must be relatively rare (upper limit of ~10%). 

CONCLUSION 

We have investigated the detection limits for sub- stellar and planetary companions to 
white dwarfs using in the WASP survey. We found that Mercury-sized bodies at small or- 
bital radii can be detected with good photometric data even in the presence of red noise. 
For smaller bodies red noise in the light-curves becomes increasingly problematic, while 
for larger orbital periods, the absence of significant numbers of in-transit points, signifi- 
cantly decreases our detection sensitivity. Application of our modified BLS algorithm to 
search for companions to WDs in our sample of 194 stars available in the WASP archive, 
did not reveal any eclipsing or transiting sub-stellar or planetary companions. We have 
used our results, to place upper limits to the frequency of sub-stellar and planetary com- 
panions to WDs. While no useful limits can be placed on the frequency of Mercury-sized 
or smaller companions, slightly stronger constraints can be placed on the frequency of 
BDs and gas giants with periods < 0.1 — 0.2days, which must certainly be relatively 
rare (< 10%). More stringent constraints would requires significantly larger WD sam- 
ples. Our key conclusion from simulations and analysis, using WASP data, suggests that 
photometric precision is of secondary importance compared to a high cadence and con- 
tinuous coverage. The short duration of eclipses and transits of WDs compared to the 
WASP observing cadence, appears to be the main factor limiting the transit detection 
rate in a survey optimised for planetary transits of main sequence stars. Future surveys 
such as Pan-STARRS and LSST will be capable of detecting tens of thousands of WDs. 
However, we emphasise the importance of high cadence and long baseline observation 
when attempting to detect the signature of close, eclipsing and transiting sub-stellar and 
planetary companions to WDs. Space missions such as COROT, Kepler (see Di Stefano 
et al. 2010) and, especially, PLATO may therefore be better suited to a survey of white 
dwarfs as they deliver uninterrupted coverage at high cadence and exquisite photometric 
precision (~ 10~ 4 — 10~ 5 ) and could at least in principle detect the transits of asteroid- 
sized bodies across a white dwarf. 
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